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The Heavy Hybrid Spectrum from NRQCD and the Born-Oppenheimer 
Approximation* 
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The spectrum of heavy-quark hybrids is studied in the leading Born-Oppenheimer (LBO) approximation and 
using leading-order NRQCD simulations with an improved gluon action on anisotropic lattices. The masses of 
four hybrid states are obtained from our simulations for lattice spacings 0.1 fm and 0.2 fm and are compared 
to the LBO predictions obtained using previously-determined glue-excited static potentials. The consistency of 
results from the two approaches reveals a compelling physical picture for heavy-quark hybrid states. 



> 

in 
O 



X 



1. Introduction 

QCD predicts the existence of hybrid states 
which contain excited gluon fields. Hybrid 
mesons with heavy bb quark pairs are the most 
amenable to theoretical treatment. They can be 
studied not only directly by numerical simulation, 
but also using the Born-Oppenheimer expansion 
which is our primary guidance for the develop- 
ment of a simple physical picture. The Born- 
Oppenheimer picture was introduced for the de- 
scription of heavy hybrid states in Refs. 0] and 
was applied using hybrid potentials first calcu- 
lated in lattice QCD in Rcf. . In this study, we 
work to leading order in the expansion and ne- 
glect higher-order terms involving spin, relativis- 
tic, and retardation effects. We test the accuracy 
of the Born-Oppenheimer approach by compari- 
son with high-precision results from simulations. 

2. NRQCD Simulation 

Iir our NRQCD simulations, we use the im- 
proved gauge-field action of Rcf. [|| on anisotropic 
lattices with as/ at = 3. The simulation parame- 
ters are listed in Table |l|. We set the aspect ratio 
using the tree level anisotropy in all our calcula- 
tions. The time symmetrized equation with sta- 
bility parameter n = 2 is used for the evolution of 
the nonrelativistic propagators. The Hamiltonian 



Table 1 

Simulation parameters for the two lattices. 





(3.0,3) 


(2.6,3) 




0.500 


0.451 


lattice 


15^x45 


10^x30 


# configs, sources 


201,16080 


355, 17040 


ro/as 


4.130(24) 


2.493(9) 


asMo 


2.6(1) 


3.9(1) 



*Talk presented by K.J. Juge. 



consists of the leading kinetic energy operator and 
two other operators to remove 0{at) and 0{a1) 
errors. 

Our meson/hybrid operators are constructed as 
follows. First the spatial links are APE-smeared. 
Then, these links are used to Jacobi-smear the 
quark (antiquark) fields and to make the (clover- 
leaf) chromomagnetic field. The specific opera- 
tors used to project out each channel are listed 
in Table ||. Note that we use 4 operators (with 
p = 0, 1, 2 and 3) in the ^ and 1~~ channels to 
get radially excited states. 

The bare quark mass is set by matching the 
ratio R = M|i„/(Eip-Eis), where M^^^ is the 
kinetic mass of the IS state, to its observed value 
21.01(6). Several low statistics runs using a range 
of quark masses were done in order to tune the 
quark mass. From the results of these runs, we 
estimate that the uncertainty in tuning the quark 



Table 2 

The meson/hybrid spin-singlet operators used in 
each channel. 



jPC 



Degeneracies 



Operator 



0-+ S 

1+- P 
Hi 

1++ H2 
0++ H, 



0,1,2++ 
0,1,2-+ 
0,1,2+- 
1+- 



T — — 



A(2) 
A 

>t B ^(2)]" 

Bx A -0 
xt B-A 



mass is about 5%. The scale is set using the av- 
erage of the IS-IP and 1S-2S sphttings. We find 
r-i = 450(15) MeV. 

The masses in the 1"' , 1++, and 0++ channels 
are extracted by fitting the single correlators to a 
single exponential for sufficiently large t. In each 
of the ^ and 1 channels, the variationally 
optimized correlator matrix is fitted to obtain the 
radially excited states. The results are shown in 
Fig. 0. 



3. Born-Oppenheimer Approximation 

In the Born-Oppenheimer approach, the hy- 
brid meson is treated analogous to a diatomic 
molecule: the slow heavy quarks correspond to 
the nuclei and the fast gluon field corresponds 
to the electrons. First, one treats the quark Q 
and antiquark Q as spatially-fixed color sources 
and determines the energy levels of the excited 
gluon field as a function of the QQ separation 
r; each of these excited energy levels defines an 
adiabatic potential ^gQ(i')- The quark motion 
is then restored by solving the Schrodinger equa- 
tion in each of these potentials. 

The energy spectrum of the excited gluon field 
in the presence of a static quark-antiquark pair 
has been determined in previous studies [|. The 
three lowest-lying levels are shown in Fig. H. 
These levels correspond to energy eigenstates of 
the excited gluon field characterized by the mag- 
nitude A of the projection of the total angular 
momentum Jg of the gluon field onto the molec- 
ular axis; by ry — ±1, the symmetry quantum 
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Figure 1. Energy splittings with respect to 
the IS state in GeV. The boxes are the Born- 
Oppenheimer spectrum where the height of the 
boxes show the effect of (5 < >= 1. 



number under the combined operations of charge 
conjugation and spatial inversion about the mid- 
point between the quark and antiquark; and for 
A = states, a refiection in a plane containing 
the molecular axis. 

Given these static potentials, the LBO spec- 
trum and wavefunctions are easily obtained by 
solving the radial Schrodinger equation with a 
centrifugal factor (Ll ) = L(L -hi) - 2A2 ^ (J^) 
where Lqq is the orbital angular momentum of 
the quark-antiquark pair and L = + J^. 

For the S+ potential, (Jg) = 0. The angular 
momentum of the excited glue is constrained by 
(Jg) > A(A + 1). For the IIu levels, we attribute 
the lowest value and for the S~ levels, we take 
(Jg) = 2 which is suggested by models of QCD. 
We show the effect of the uncertainty in this quan- 
tity by changing the value by one unit. The par- 
ity P and charge conjugation C of each meson is 
given in terms of L and S by P = e (_i)L+a+i 
and C = e 77 (-1)'"+'^+^, where L > A and e = 1 
for E+, e = -1 for E", and e = ±1 for A > 0. 

Results for the LBO spectrum (splittings) of 
conventional bb and hybrid bgb states are shown 



Table 3 

Comparison between LBO and NRQCD split- 
tings in GeV. The first error in LBO is due to 
a 5% change in mt and the second error is an 
estimate of the uncertainty due to ^ < >. 







LBO 


NRQCD 


% diff 


IP 


- IS 


0.392(1) 


0.432(4) 


9(1) 


2S 


" IS 


0.553(5) 


0.586(5) 


6(1) 


3S 


- IS 


0.957(9) 


0.950(8) 


1(1) 


Hi 


- IS 


1.421(8)(37) 


1.479(24) 


4(3) 


H2 


- IS 


L421(8)(37) 


1.517(58) 


6(5) 


H.3 


- IS 


L711(5)(63) 


1.808(25) 


5(4) 


H'l 


- IS 


L697(2)(36) 


1.892(30) 


10(2) 



in Fig. |l|. The scale vq is taken from the NRQCD 
simulation. The heavy quark mass Mb is tuned to 
reproduce the experimentally-known T(1S) mass: 
Mt — 2Mb + Eq, where Eg is the energy of the 
lowest-lying state in the potential. Level split- 
tings are insensitive to small changes in the heavy 
quark mass. For example, a 5% change in Mb re- 
sults in changes to the splittings (with respect to 
the IS state) ranging from 0.1 — 0.8%. 

The applicability of the leading Born Oppen- 
heimer approximation relies on the smallness of 
the retardation effects. We can determine the size 
of this effect by comparing the two approaches 
taken here since the NRQCD Hamiltonian only 
differs from the LBO Hamiltonian by the inclu- 
sion of the p • A coupling between the color charge 
in motion and the gluon field. The difference in 
splittings are tabulated in Table ^ and shown in 
Fig. 1^. In the comparison, the scale ambiguity 
has been removed and the lattice artifacts were 
found to be small (< 4%). Furthermore, the 
splittings are insensitive to small changes in the 
&-quark mass. We thus conclude that the small 
differences (< 10%) from the two approaches is 
in fact due to retardation, hence validating the 
Born-Oppenheimer expansion. 

4. Conclusion 

The masses of 4 hybrid states were determined 
in leading order NRQCD in the quenched approx- 
imation. The magnitude of the retardation ef- 
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Figure 2. Wavcfunctions and potentials for the 
various hybrid/meson states. 

fects were found to be small validating the Born- 
Oppenheimer expansion. 

The level splittings in the LBO approximation 
were found to be rather insensitive to the b quark 
mass. The LBO wavcfunctions revealed the spa- 
tial largeness of the hybrid meson states and also 
indicated that the finite volume effects in our sim- 
ulations should be negligible. Note that our LBO 
spectrum of hybrid mesons based on adiabatic 
surfaces calculated in the quenched approxima- 
tion will very likely differ from the true spectrum 
due to our neglect of sea quark effects. The in- 
clusion of such effects remains an important chal- 
lenge for the future. 
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